Generalized Bounds on Majoron-neutrino couplings 



R. Tomas, H. Pas Q J. W. F. Valle 

Institut de Fisica Corpuscular - C.S.I.C., Departament de Fisica Teorica - Univ. de Valencia 
Edifici Instituts d'Investigacio - Apartat de Correus 2085 - 46071 Valencia, Spain 



Abstract 

O 

o 

CN ' We discuss limits on neutrino-Majoron couplings both from laboratory experi- 

ments as well as from astrophysics. They apply to the simplest class of Majoron 
1^ . models which covers a variety of possibilities where neutrinos acquire mass either via 

■ a seesaw-type scheme or via radiative corrections. By adopting a general framework 

CN I including CP phases we generalize bounds obtained previously. The combination 

. of complementary bounds enables us to obtain a highly non-trivial exclusion region 

O ' ill the parameter space. We find that the future double beta project GENIUS, to- 

, gether with constraints based on supernova energy release arguments, could restrict 

neutrino-Majoron couplings down to the 10"'' level. 

^ : 1 Introduction 



X 



The confirmation of the zenith-angle-dependent atmospheric neutrino deficit by the Su- 
perkamiokande experiment generally has been understood as first significant hint for neu- 
trino masses and thus particle physics beyond the standard model The other long- 
standing puzzle of particle physics is the deficit of solar neutrinos 0. Altogether they 
constitute the most important milestone in the search for phenomena beyond the Stan- 
dard Model (SM), indicating the need for oscillations involving all three active neutrino 
species P]. The mounting experimental activity in this field promises a bright future for 
neutrino physics which may prove to be a most valuable source of information on the struc- 
ture of a more complete theory underlying the standard model of particle physics. 

An elegant way to introduce neutrino masses is via the spontaneous breaking of an 
ungauged lepton number symmetry through a non-zero SU{2) U{1) singlet vacuum ex- 
pectation value (VEV) of a scalar field. This may be implemented in conventional |Q, ^ as 
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well as supersymmetric models 0. The couplings of the corresponding Goldstone boson, 
generically called Majoron and denoted by J, are rather model-dependent [0. Here we con- 
sider the simplest class of Majoron models, where the Majoron- neutrino coupling matrix 
gf- oc TJiij is proportional to the neutrino mass matrix p|, so that in the mass eigenstate 
basis diag (mi,m2,'m3) the Majoron neutrino couplings Qi are diagonal, to lowest order 
approximation [], 

gfl ^ h9r- (1) 

This covers a variety of possibilities including both seesaw-type as well as radiative 
models [Q. 

Limits on this quantity obtained from laboratory experiments searching for Majoron- 
emitting pion or kaon decays are rather weak, with the exception of double beta decay [jll] . 



On the other hand Majoron emitting neutrino decays affect the expected neutrino luminosity 
and spectra which are constrained by the observed signal from SN1987A, providing stringent 
restrictions . While the limits of laboratory experiments on rare decays are given in the 



weak basis, bounds from processes in supernovae occur in a dense medium and are expressed 
in the medium eigenstates (see below). In the present work we discuss the correlations of 
the different limits and their translation into the mass basis, extending the earlier paper 
ref. |T^. In the next section we derive the expressions for medium and weak eigenstates, 
following |jl2|, |l^. In section 3 we review the bounds obtained from the supernova SN1987A 



using various considerations. In contrast to ref. [|12| here we include the study of the effects 
associated to the Majorana CP violating phases present in theories of massive neutrinos |n, 
p!5| . Moreover we investigate (Section 4) the recent bounds from neutrinoless double beta 
decay as well as those that could be attained in future experiments such as GENIUS []TB 
The resulting exclusion plots are discussed in section 5 in the mass basis. 



2 Neutrino mixing in three bases 

For neutrinos propagating through a medium one has to deal with three kinds of eigen- 
states: Flavor eigenstates Va, mass eigenstates vf'''^ with masses rrij, and, depending on the 
environment, medium eigenstates i'f^''\ The flavor eigenstates are defined as 

i^a = Ef^-^f"^ (2) 

i 

^This proportionality may be avoided in more complex models, such as those in ref. |9[ 
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and the medium eigenstates are 

'"'^"'^-EM'^'^- (3) 



Here the superscript {hi) = ±1 refers to the hehcity of the state. In the general case it is 
impossible to diagonalize simultaneously the mass and potential terms. Thus one has to 
solve the field equations in detail. In a two-component field formalism, where a left-handed 
four-component field u expressed in the chiral representation of the 7 matrices fl^ is related 



to the corresponding two-component field by z/J = (0^, 0) [|T3] | , the Lagrangian can be 
written in the mass basis as 

m. 



^E^?.?('/'f^20,+ 01^20*) (4) 



where (J2 and a denote Pauli matrices. Here the free Langrangian Cq and £med describe 
the propagation in vacuo and the effects of matter described by the potential matrix Vij, 
respectively, whereas £int takes into account the presence of neutrino-Majoron interactions 
which may lead to decays. One has now to consider the decays ^iXpi) — > ^^\Pj) + J{(i)i 
where v'lXpi) and v^\pj) are energy-eigenstate Majorana neutrinos that propagate in matter 
with four- momenta Pi = {E^',pi) and pj = {Ej\pj), and helicity hi and hf respectively. 
In order to obtain these energy-eigenstates one has to take £0 + £med and calculate the 
resulting field equations. 

{ido -iax V)0i(x) + miia2(j)*i - E = 0- (5) 

One solves these field equations by expanding the fields 0i(x) as superpositions of plane- 
wave spinors with definite helicity, |jT3|, IT] 



where a{p) and j3{k) are helicity eigenstates and and A^* denote positive and nega- 
tive frequency components of the field under consideration. One should now substitute 



^The notation here coincides to the one of Q up to a factor of ; 
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Table 1: Medium eigenstates vf, equivalent weak eigenstates in the limit \Vaa\ ^ "^f/(2p) 
and their potential energy. The rotation in the — subspace is parametrized by C23 = 
cos6'23 and S23 = sin ^23, where the arbitrary argument has been chosen to coincide with 
the mixing angle ^23 in vacuo (from [0]). 



this expression in the equations whose diagonalization would give rise to the desired 
eigenstates. It can be shown, though, that for relativistic neutrinos the positive-frequency 
components decouple from the negative-frequency ones and the energy eigenstates obtained 
in this way result to be the same as those obtained from the diagonalization of the usual 



MSW equation |jT3[, which can be stated as 

Here Hlf (p + mf / {2p))5ij and Vap is the potential matrix in the weak basis, 

/ V^^V.r \ 



(7) 



a/3 



Vc + Vn 





Vn 

v^^ ) 



(8) 



The potentials induced by the charged and neutral currents are Vc = \/2hGFnBiYe + 
and Vat = \/2hGFnB (~|^Af + ^i^e); where = {rii — n-j)/nB and ub is the baryon density. 
Diagonalizing H"^^^ + UVU^ yields the medium eigenstates 



~{h) 



In the three-flavor neutrino case the mixing matrix U can be parametrized as U 



U23U13U12U0, where the matrices Uij = t/, 



angle 9ij and Uq includes possible CP-violation effects \15 



perform the rotation in the ij-plane by the 
In the following we will 



assumme ^13 = 0, motivated both by detailed fits of the present solar and atmospheric 
neutrino anomalies as well as by the reactor results of the Chooz experiment ||18[. This 
simplifies the mixing matrix to Ua = Uai^i = U23Ui2Uoi^i and allows us to set = 9q 
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and 623 = Oatm- Now for light neutrinos near the neutrinospheres in supernovae the condition 
\Vaa\ ^ "^f/(2p) holds and, since in the weak basis the potential is diagonal, the medium 
states can be identified with the weak ones up to an arbitrary rotation in the z/^ — z/^ 
subspace. In order to simplify the expressions we exploit this freedom by choosing this 
arbitrary rotation angle to coincide with —^23? see Table |l|. This allows us to identify the 
coupling matrix in the medium basis with the one in the weak basis up to the rotation 

Taking now into account that g^'^ = U'^ U and substituting the explicit expressions 
for the Uai matrices relating mass and weak eigenstates, one gets the following expression, 
~9ij = U12U* gf^ UlU^2^ or explicitly 

' 



g 



gee 
gii'e 



113 
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This choice of the rotation in the z/^ — z/,- subspace leads to a relation between medium and 
mass eigenstates characterized only by the solar angle 6q and by the Majorana CP violating 
phase 6 [H, 15]. Using the definitions Am^2 = and Am23 = Am^j^ together with 

the assumptions in ^ one can easily translate the bounds obtained in the weak or medium 
basis into the mass basis and, in addition, express them in terms of only two independent 
parameters, for instance (mi, (71) via 



g2 = gi 



\ 



1 



Ami 

2 

mi 



5^3 = gi 



1 



Am| + Am^ 



atm 



2 

mi 



fill 



.(10) 



3 Supernova bounds 



There is a variety of different arguments based on supernova physics which lead to restric- 
tions on neutrino properties. Processes involving Majoron-neutrino couplings may prevent 
a successful explosion as well as substantially affect the observed neutrino spectra. A cru- 
cial feature to notice is that the effective mass induced by the interactions of neutrinos 
with background matter breaks the proportionality between the neutrino mass matrix and 
the neutrino-Majoron coupling matrix . This follows from the fact that the thermal 
background in the supernova environment consists only of particles of the 1st generation, 
thus distinguishing the electron fiavour from the others. We now describe three different 
arguments used |]T2[ in order to restrict the relevant parameters. 
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3.1 Constraints from Neutrino spectra 

The idea behind this bound is that Majoron-induced transitions between the neutrino 
flavors could change the energy spectra of the single flavors. At the typical temperatures 
of the SN core u^ .^ only interact with the medium via neutral currents giving rise to a 
smaller cross section than that corresponding to the electron (anti)neutrinos, which feel 
both neutral and charged currents. Since the opacity of the heavier i/^ flavors is smaller 
than for the z/g, their energy-exchanging reactions freeze out in the denser region of the 
protoneutron star, leading to lower spectral temperatures of z/g compared to z/^^^. This 
expected spectrum can be distorted due to the decays vj + J. Besides the effects of 

such decays one has to keep in mind the possible oscillation which neutrinos could undergo 
along their journey to the Earth. In order to consider both aspects we have defined the 
effective survival probability as 

N = Ndecay X A^..^, (12) 

where Ndecay stands for the survival probability of a emitted from its energy sphere and 
can be computed as 

Ndecayii>t) ^ exp | - H dr ^ r(z>f ^ Z>J + J) i . (13) 

Within our relativistic approximation, the helicity-flipping neutrino decays rate are given 
by 

nDt^yJ + J)=^-^{V,-Vj). (14) 

Coming to the oscillation term, the corresponding neutrino survival probability Noso will de- 
pend on the neutrino mixing angle and squared mass difference. We will analyse separately 
the different solutions of the solar neutrino problem namely small-angle MSW (SMA-MSW), 
large-angle MSW (LMA-MSW), LOW-MSW and the just-so case. Details about the present 
status and required parameters of the various solutions can be found in the global analysis 
of neutrino data presented in ref. |p[. Such a study favours a rather small value for the 



angle ^13, mainly because of data from reactors In the first three cases neutrinos will 
propagate through the supernova environment adiabatically. Therefore they will emerge as 
energy eigenstates, which in vacuum coincide with the mass eigenstates, without any oscil- 
lation occuring on the way from the SN to Earth. If one takes into account that neutrinos 
have to traverse a distance, rf, of matter in the Earth to reach the detectors, Kamiokande 
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and 1MB, one gets the following expression for their survival probability ^ 

Nose = 1 - {sin^ 9q - sin 29^ sin(2^Q - 29 J sin^ ( ^ ) }, (15) 

where Im and 9^ denote the oscillation length and the mixing angle in matter, respectively. 
As has been previously noted for the simplest choice ^13 = one has that, besides the fact 
that and u^- behave the same way in the supernova, the conversion v^i will be 

the only oscillation involving electron (anti)neutrinos, allowing us to set the angle which 
characterizes their mixing to 6'©. 

In the vacuum solution case the neutrinos emerge essentially as flavor eigenstates which 
then oscillate on their way to Earth. Therefore one has 

iVosc = l-^sin2 2^0. (16) 

In order to get information on the coupling constants we will conservatively require that at 
least half of the initial electron antineutrinos emerging from the SN1987A survive, > 0.5, 
accounting for the rough agreement between the expected and the detected SN1987A signals. 
In order to analyse the implications of this restriction one must generalize the simplest 
argument used in pOf since neutrinos may loose energy as a result of majoron decays. 

This allows us to get some limits on the coupling parameter of the order of (71 ((72) < 
few X 10~^ from the first three solutions to the solar neutrino problem. For the case of 
vacuum oscillations, though, the solution is already disfavored by the SN1987A data even 
in the absence of neutrino decays pO|. Though they may narrow it down considerably, the 



above arguments do not totally close the allowed window of neutrino-Majoron couplings, 
neither for the SMA, LMA nor LOW solutions, even for a supernova in our milky way. 

3.2 Constraints from Majoron luminosity 

This bound is based on the observation that neutrino decays into Majorons could supress 
the energy release contained in the neutrino signal. Under the assumption of small Ve — 
mixing the neutrino signal observed in SN 1987A is in good agreement with numerical 
computations of the total binding energy released in a supernova explosion. An analysis of 
decay and scattering processes involved yields the exclusion region [|r^ 

3 X 10"^ < \~g,j\ < 2 X 10"^ (17) 

For \gij\ values smaller than 3 x 10^^ the Majoron neutrino coupling becomes too small to 
induce any effect. On the other hand for \gij\ > 2 x 10^^ Majorons get trapped in the core 
and do not contribute to the energy release. 
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Another point to observe is that CP violating phases affect these hmits. This follows 
from the appearance of the phase 6 in explicit form of the Majoron neutrino coupling 



constants given in eq. ([10|) . In order to eliminate such an explict CP phase dependence 
when translating the limit on \gij\ into the mass basis we have analyzed for each term 
of eq. (|10D the excluded region for different values of 6 and subsequently considered the 
intersection of the resulting excluded regions. This conservative procedure allows us to rule 
out part of the parameter space irrespective of the value of the CP phase. 

As an example we illustrate in Figs. 1 and 2 the regions excluded for the LMA MSW 
solution to the solar neutrino problem. The luminosity bound can be described in two 
steps. In the first one we take one Qij from eq. (|10[) and by means of eq. ([TI|) we write it 
in terms of gi — mi. This way we obtain an expression for the energy loss which depends 
explicitly upon the CP phase 6. Now, by varying that CP phase the bound given in eq. 
(0) is translated into different ruled out regions. We show in Fig. 1 the resulting bound 
on \gee\ assuming two extreme cases, 6 = (solid lines), and 6 = it /2 (dotted lines). Notice 
that for the latter case the bound disappears because of a cancellation between the two 
terms in \gee\- In order to remove the phase we therefore consider the intersection as the 
most conservative choice. 

Now turn to the implications of the luminosity bound to the other components of 
the Majoron neutrino coupling matrix elements. Once we have obtained those intersecting 
regions for each cjij we simply take the union of them, giving rise to a final highly non-trivial 
exclusion region, as can be seen in Fig. 2. 

It is important to notice that the shape of such regions is characterized by the values of 
the square root of Am| and Arn^j^. Let us first consider gij with i,j ^ 3. In this case only 
Amg appears in eq. ([11]), so that for mi ^ Am^ one has g^ ~ g\ giving rise to a vertical 
line with no dependence on mi, as noted in the figures. In contrast, for mi <^ Am^ one 



_ "0 
has ~ mi 'with a explicit dependence on mi which strengthens the bound for lower 



mi values. Let us now consider the limit coming from g-^. In this case the characteristic 
mass scale is always given by Am^^^, eq. {^), irrespective of the particular solutions to 
solar neutrino problem that one may wish to consider. As a result, for the LOW and 
(quasi)-vacuum cases the difference between Am^^^ and the solar mass scale AitLq is so 
large that two branches appear. This explains the two branches observed in figures 5 and 
6 corresponding to the LOW and (quasi) -vacuum solutions, respectively (section 4). 

Concerning the SMA solution in fig. 4 the main changes arise from the bounds on \gee\ 
and \geii'\- In the expression of gee in eq. (10) the contributions of gi and g2 may cancel for 
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Figure 1: Excluded regions from the Majoron luminosity requirement 3 x 10"'' < \gee\ < 
2 X 10"^, for two extreme CP cases, 6 = (solid lines) and 5 = n/2 (dotted lines), in the 
nil — gi plane. LMA parameters sin^(2^) = 0.6 , Am^ = 1 x 10~^ eV^ are assumed. 



a phase 5 ~ 7r/2, see fig. 1, as long as these contributions are of comparable magnitude. To 
fulfill this requirement, a smaller admixture of g2 in g^e, as happens for the SMA solution, 
has to be compensated by a larger ratio of g2/gi, corresponding to smaller masses mi. This 
is producing a small hollow in the bound at mi ~ 5 x 10~®. 

The bound on \ge^_i'\ is responsible for the sharp peak at the right edge of the excluded 
region, as can be seen in fig. 2, where all Majoron luminosity bounds are shown explicitly. 
Here the conservative upper bound on \gefi'\ is obtained, when 6 = 0, corresponding to 
cancellation of the gi and g2 contributions in the large mi asymptotics, see eq. (10). 
Correspondingly the right border of the excluded region is obtained for 6 = 7r/2, where 
\geti'\ becomes constant for large mi. The intersection gives rise to the peak. For smaller 
values of the mixing, as obtained for the SMA solution, the expression for \gefi'\ in eq. (10) 
is fulfilled for corresponding larger values of gi. This shifts the exclusion region to the right, 
making it more visible in fig. 4. 

Before concluding this section we mention the constraints on Majoron- neutrino coupling 
parameters which arise from the collapsing phase. The idea behind this bound is that a 
change in the trapped electron fraction could prevent a successful explosion process. At 
the end of their life massive stars become unstable and, when the iron core reaches the 
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Figure 2: Excluded regions independent of CP phase from the Majoron luminosity require- 
ment on |(7ee|, (solid line), \geij.'\ (dashed line), Ig^'n'l (dotted lines) and {Qt't'I (dash-dotted 
line), in the mi — gi plane. LMA parameters sin^(20) = 0.6 , AttIq = 1 x 10~^ eV^ are 
assumed. 



Chandrasekar limit, they implode. Once the nuclear density is reached, a shock wave forms 
at the core and propagates outwards, turning the implosion into an explosion. The strength 
and propagation of this shock is sensitive to the trapped electron fraction Y^^ = + Y^,^, 
which can be erased by neutrino decays ^ Vf, + J. Requiring Yi(tf,ounce) > 0.375 leads 



to a limit of [0 

\gee\ < 2 X 10"^ (18) 

However to the extent that current supernova models do not fully account for the explosion 
mechanism this limit should be taken only as indicative for the moment. 



4 Neutrinoless double beta decay 

The only laboratory experiment, which is competitive with the supernova bounds, is neu- 
trinoless double beta decay. This decay corresponds to two single beta decays occuring in 
one nucleus and converts a nucleus (Z,A) into a nucleus (Z+2,A). Limits on the Majoron 
emitting mode 

iX-^i^,X + 2e- + ^ (19) 
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are given by two types of experiments. In geochemical experiments the half-life limit is 
derived from relative abundances of nuclear isotopes found in the earth [B3| : 



\9e 



< 3 X lo-^ 



(20) 



However, half-life determinations vary by more than a factor of three. 

The best direct laboratory limit (less stringent but more reliable) from the Heidelberg- 



Moscow experiment |23] is based on a likelyhood fit to the continuous electron spectrum: 

\gee\ < 8 X 10-^ (21) 



Future projects such as GENIUS |jT^ and EXO p4| aim at considerable improvements in 
the sensitivity. A very rough estimation of the sensitivity of GENIUS It is based on the 



background simulation in [0, where a background improvement in the interesting energy 
range of a factor ~ 1000 has been obtained. Since in the Heidelberg-Moscow experiment 
the Majoron-neutrino coupling bounds are dominated by the systematical error of the back- 
ground simulation, a considerable reduction of the background by a factor of B will reduce 
the limit on the Majoron-emitting double beta decay half life by ~ ^/B and the coupling 
constant limit accordingly by r^i/B. This implies a reach of sensitivity down to \gee\ ~ 10^^, 
which could bridge most of the gap existing between the more reliable limits derived from 
super novae. 



5 Discussion and conclusions 

In figures 3 to 6 we present the limits on Majoron-neutrino couplings in terms of mi — gi 
corresponding to the various solutions of the solar neutrino problem. In Fig. 7 we display 
the results for the LMA solution also in terms of the equivalent m2 — g2 variables. This 
representation has been selected for convenience and generality. By further specifying the 
underlying model for lepton number violation one can re-express our results in terms of the 
lepton number breaking VEV, which will provide also useful information for model-builders. 

Regions which are excluded by supernova arguments are denoted by the rhombical 
pattern (obtained from Majoron luminosity) and by the vertical lines (obtained from the 
neutrino spectra). Also shown are the regions excluded from neutrinoless double beta decay 
(horizontal lines). The excluded region from Majoron luminosity is a superposition of the 
bounds on gn, gi2, g22, and ^33, where always the most conservative limits for various CP 
Majorana phases have been used. Due to the expressions for the helicity-fiipping neutrino 
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decays the bound obtained from the neutrino spectra turns out to be independent of the 
CP phase. For neutrinoless double beta decay one can have a cancellation of the coupling 
constants gi and g2- The expected sensitivity of the GENIUS experiment is shown as a 
dashed line. It is easy to see that GENIUS could be able to bridge almost the whole gap 
between the different supernova constraints. Also an upper bound mi < 2.3 eV from tritium 
beta decay is displayed. 

The limits obtained in this paper apply to the simplest class of models where neutrino 
masses arise from the spontaneous violation of lepton number. Such Majoron models which 
cover a wide and attractive class including both models where the smallness of neutrino 
masses follow from a seesaw scheme, as well as those where it arises from the radiative 
corrections. 

Both neutrinoless double beta decay as well as supernova physics arguments provide 
stringent limits on Majoron-neutrino interactions. In the present work we have discussed 
these limits and their translation into the mass basis. Generalizing previous papers |12| we 



have now taken into account the effect of CP violating phases, which play a crucial role in 
the neutrinoless double beta decay limits. Depending on the solution of the solar neutrino 
problem and the absolute mass scale in the neutrino sector the constraint from the supernova 
energy release (Majoron luminosity argument) exludes Majoron-neutrino couplings in the 
wide range of 10^^ — 10~^. Upper bounds have been obtained from neutrinoless double 
beta decay and the SN87A neutrino spectra. An estimate of the potential of the future 
double beta projects such as GENIUS suggests the possibility to bridge almost the whole 
gap separating the excluded areas and either to establish Majorons with couphngs around 
a few 10^^ or to restrict neutrino-Majoron couplings down to 10^''. 

Last, but not least, let us mention that the propagation of neutrinos produced in the 
solar interior follows essentially the MSW picture, while any possible effect of decays would 
happen in vacuo through a non-diagonal neutrino- majoron coupling which is absent in the 
simplest models considered here [Q. Even in more complex models |l^ where such non- 
diagonal neutrino-majoron couplings exist in vacuo, one can see that for such small values 
of the neutrino-majoron coupling strengths indicated by supernova and neutrinoless double 
beta decay, it is rather unlikely that they can play any role whatsoever in the solar neutrino 
problem . 
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Figure 3: Bounds on Majoron models in the gi — mi plane for the case of the LMA solution 
(Am| = 10"^ eV^ and sin^ ~ 0.6). 
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Figure 6: Bounds on Majoron models in the gi — mi plane for the vacuum solution of the 
solar neutrino problem {AttIq = 10~^° eV^ and sin^ 29q ~ 0.9. The bounds from neutrino 
spectra don't apply, since already pure neutrino oscillations implied by the vacuum solution 
parameters lead to a contradiction with the observed neutrino spectrum of SN1987A. 
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Figure 7: Bounds on Majoron models expressed in terms ofthe g2 — m2 for the case for the 
LMA solution of the solar neutrino problem. 

Acknowledgement 

We thank M. Hirsch for useful discussions. This work was supported by Spanish DGICYT 
under grant PB98-0693, by the European Commission TMR contract HPRN-CT-2000- 
00148 and by the European Science Foundation network grant N. 86. H. P. was supported 
by TMR contract ERBFMRX-CT96-0090 and DOE grant no. DE-FG05-85ER40226 and 
R.T. by a grant from the Generalitat Valenciana. 



15 



References 



[1] H. Sobel, \http://nu2UUU. sno. Laurentian. c(\ T. Toshito, fittp://www. ichep2UUU. rL. ac. uk . 
Y. Fukuda et al. [SuperKamiokande Collaboration], Phys. Rev. Lett. 81 (1998) 1562; 
Phys. Rev. Lett. 82 (1999) 2644 ||hep-ex/9812"0Tl . Y. Fukuda et al, Phys. Lett. B433, 
9 (1998); Phys. Lett. B436, 33 (1998); For a recent analysis see: N. Fornengo, M.C. 
Gonzalez-Garcia, J. W. F. Valle, Nucl. Phys. B580 (2000) 58-82 ||hep-ph/000:ZTTr | 



[2] Y. Suzuki, talk at XIX International Conference on Neutrino Physics and Astrophysics, 
Sudbury, Canada, \http://nu2000. sno. laurentian. cat T. Takeuchi, talk at the XXXth In- 
ternational Conference on High Energy Physics, ICHEP 2000 http:/ / www. ichep2000. - 
rL.ac.ul\ . Homestake Collaboration, B.T. Cleveland et al., Astrophys. J. 496, 505 
(1998); R. Davis, Prog. Part. Nucl. Phys. 32, 13 (1994); K. Lande, pttp://nu2UUU.sno. 



laurentian. cc^ , SAGE Collaboration, J.N. Abdurashitov et al., Phys. Rev. C60, 055801 
(1999); V. Gavrin, fittp://nu2000. sno. laurentian. c(\ . GALLEX Collaboration, W. Ham- 



pel et al, Phys. Lett. B447, 127 (1999). E. Belloti, }[ittp://nu2UUU. sno. laurentian. cc\ . 
For a recent analysis see: M. C. Gonzalez-Garcia, P. C. de Holanda, C. Pena-Garay 
and J. W. F. Valle, Nucl. Phys. B573 (2000) 3 ||hep-ph/9906469 | 



[3] M. C. Gonzalez-Garcia, M. Maltoni, C. Pena-Garay and J. W. F. Valle, Phys. Rev. D 
63 (2001) 033005 ||hep-ph/0009350[ . 



[4] Y. Chikashige, R.N. Mohapatra and R.D. Peccei, Phys. Lett. B 98 (1981) 26 

[5] Two basic types of radiative neutrino mass models were considered in A. Zee, Phys. 
Lett. B 93, 3890 (1980) and K. S. Babu, Phys. Lett. B 203, 132 (1988). For Majoron 
embeddings of these models see, for example, A. S. Joshipura and J. W. F. Valle, Nucl. 
Phys. B 397, 105 (1993). 

[6] A. Masiero and J. W. F. Valle, Phys. Lett. B251 (1990) 273 

[7] J.W.F. Valle, Gauge Theories and the Physics of Neutrino Mass, Progr. Part. Nucl. 
Phys. 26 (1991) 91-171 

[8] J. Schechter and J. W. F. Valle, Phys. Rev. D 25 (1982) 774. 

[9] M. C. Gonzalez-Garcia and J. W. F. Valle, Phys. Lett. B 216 (1989) 360. 

[10] G. B. Gelmini and J. W. F. Valle, Phys. Lett. B 142 (1984) 181. 



16 



[11] M. Hirsch, H.V. Klapdor-Kleingrothaus, S.G. Kovalenko, H. Pas, Phys. Lett. B 372 
(1996) 8; M. Hirsch. H. V. Klapdor-Kleingrothaus, B. Maier, H. Pas, in Proc. Int. 
Workshop on Double Beta Decay and Related Topics, Trento, 24.4.-5.5.95, World Sci- 
entific Singapore, Ed.: H.V. Klapdor-Kleingrothaus and S. Stoica 

[12] M. Kachelriess, R. Tomas, J.W.F. Valle, Phys. Rev. D 62 (2000) 023004 ||hep-ph 



000103911 



[13] C. Giunti, C.W. Kim, U.W. Lee, W.P. Lam, Phys. Rev. D 45 (1992) 1557 
[14] J. Schechter and J.W.F. Valle, Phys. Rev. D22 (1980) 2227. 
[15] J. Schechter and J. W. F. Valle, Phys. Rev. D 23 (1981) 1666. 

[16] H.V. Klapdor-Kleingrothaus, Proc. Beyond the Desert '97; J. Hellmig, H.V. Klapdor- 
Kleingrothaus, Z. Phys. A 359 (1997) 351; H.V. Klapdor-Kleingrothaus, M. Hirsch, Z. 
Phys. A 359 (1997) 361; H.V. Klapdor-Kleingrothaus, J. Hellmig, M. Hirsch, J. Phys. 
G 24 (1998) 483; H.V. Klapdor-Kleingrothaus, L. Baudis, G. Heusser, B. Majorovits, 
H. Pas, |hep-ph/9910205| 

[17] P.D. Mannheim, Phys. Rev. D37, 1935 (1988). 

[18] M. Apollonio et al., Phys.Lett. B466 (1999) 415, |hep-ex/9907037| ; F. Boehm et al., 
|hep-ex/9gT^ 

[19] This simplified form of the lepton mixing matrix was first given in J. Schechter and 
J.W.F. Valle, Phys. Rev. D21, 309 (1980). 

[20] A.Yu. Smirnov, D.N. Spergel and J.N. Bahcall, Phys. Rev. D49 (1994) 1389. 

[21] B. Jegerlehner, F: Neubig and G. Raffelt, Phys. Rev. D54 (1996) 1194. 

[22] T. Bernatowicz et al., Phys. Rev. Lett. 69 (1992) 2341 

[23] H.V. Klapdor-Kleingrothaus et al. (Heidelberg-Moscow CoUab) [hep-ph/0103062| ; L. 
Baudis, priv. comm. 

[24] M. Danilov et al., Phys.Lett. B480 (2000) 12-18 

[25] J. N. Bahcall, S. T. Petcov, S. Toshev and J. W. Valle, Phys. Lett. B 181 (1986) 369. 
For a recent paper see, e. g. A. Bandyopadhyay, S. Choubey and S. Goswami, Phys. 
Rev. D 63 (2001) 113019 



17 



